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Abstract 

In contrast to scalar and tensor modes, vector modes of linear perturbations 
around an expanding Friedmann-Robertson-Walker universe decay. This makes 
them largely irrelevant for late time cosmology, assuming that all modes started 
out at a similar magnitude at some early stage. By now, however, bouncing models 
are frequently considered which exhibit a collapsing phase. Before this phase reaches 
a minimum size and re-expands, vector modes grow. Such modes are thus relevant 
for the bounce and may even signal the breakdown of perturbation theory if the 
growth is too strong. Here, a gauge invariant formulation of vector mode perturba- 
tions in Hamiltonian cosmology is presented. This lays out a framework for studying 
possible canonical quantum gravity effects, such as those of loop quantum gravity, 
at an effective level. As an explicit example, typical quantum corrections, namely 
those coming from inverse densitized triad components and holonomies, are shown to 
increase the growth rate of vector perturbations in the contracting phase, but only 
slightly. Effects at the bounce of the background geometry can, however, be much 
stronger. 



1 Introduction 



Linear perturbations on isotropic expanding or contracting geometries can be split in dif- 
ferent types according to their transformation properties under spatial rotations: scalar, 
vector and tensor modes. In general, this presents a convenient decomposition of general 
perturbations into different classes. But it acquires a more important role for the linearized 
Einstein's equation where the different modes decouple. Moreover, linearized gauge trans- 
formations corresponding to changes of space-time coordinates do not mix the modes which 
can thus be analyzed separately. In this article we focus on the vector mode and possible 
quantum effects in its evolution. 

Vector mode perturbations in classical cosmology decay in an expanding universe, and 
dynamically become of less interest as the universe continues its expansion. Accordingly, 
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vector modes are often ignored. This feature holds true for expanding cosmologies which 
start from a big bang or emerge from a quantum state if one assumes all modes to be of 
comparable initial magnitude. However, this assumption has to be justified. One possibility 
is to use bouncing models where the pre-history before the big bang is described by a 
collapsing phase. Then, vector perturbations can lead to significant problems because 
gauge invariant measures of vector perturbations grow. Their current size relative to that 
of scalar and tensor modes then depends on where equal sizes or other initial conditions 
are assumed. Moreover, vector modes are generated by higher order perturbations and 
subsequently grow in a contracting phase p. They can thus not be ignored altogether. 

There are several models such as string inspired pre-big bang scenarios [21 [3] or cyclic 
and ekpyrotic models [H [5] , or some models of loop quantum cosmology [6l [7] which exhibit 
a bounce. Such scenarios for an avoidance of the big bang singularity are developed mainly 
based on homogeneous models of cosmology. On the other hand, the growth of vector 
perturbations in the contracting phase indicates a possible violation of the homogeneity 
assumption when the bounce is approached as indicated by the breakdown of classical 
perturbation theory. Thus, the growth of vector perturbations not only raises questions 
regarding the validity of an homogeneity assumption but may even question the phenomena 
of the bounce itself. Q 

In this paper, we derive the vector mode dynamics in the context of cosmological mod- 
els based on loop quantum gravity [HI [121 IE] and cosmology [H] . Firstly, we present a 
systematic derivation of classical vector perturbation equations using a canonical formu- 
lation in Ashtekar variables [151 IIHI- We compute the gauge transformations property for 
the vector perturbations and then construct the corresponding gauge invariant variable. 
All this is done in a purely canonical way to outline the general procedure followed also in 
the presence of quantum corrections. 

In the following section we study possible effects of quantum corrections expected from 
loop quantum gravity. We compute requirements on quantum correction functions from 
anomaly cancellation in the quantum corrected constraint algebra. These expressions in- 
clude inhomogeneities in such a way that all symmetries of a Friedmann-Robert son- Walker 
background are broken. Thus, we are not dealing with a mini-superspace quantization even 
though inhomogeneities are restricted to the perturbative vector mode. How individual 
terms of effective constraints are related to operators in the full theory is described in [17] . 
We study the effects of quantum correction arising for inverse densitized triads in detail, 
and perform corresponding calculations for corrections from holonomies in a subsequent 
section. While current methods are not sufficient to compute full correction functions for 
all gauges, it turns out that the remaining freedom is constrained by requiring the absence 
of anomalies as a consistency condition. This provides evidence, at the effective level em- 
ployed here, that anomaly cancellation will restrict possible quantization choices of the full 

^Based on matching conditions, an evolution of vector modes through a bounce has been studied, 
e.g., in [51[S], assuming certain forms of a non-singular bouncing background geometry. This does not 
address the validity of a perturbativity assumption. Vector modes have occasionally been considered in 
inflationary scenarios such as in [10 . Since this happens in an expanding phase, vector modes decay and 
do not challenge the perturbativity assumption either. 



2 



theory. (See [T8l[T9] for analogous statements based on different principles.) Moreover, we 
demonstrate explicitly that anomaly-free set of constraints, and thus a covariant effective 
space-time picture, is possible even in the presence of non-trivial quantum corrections. As 
one application, we show that in perturbative regimes (not close to a bounce) quantum 
corrections make the growth rate of the vector mode in a contracting universe slightly 
stronger compared to the classical behavior. 

2 Canonical formulation 

In this paper, we study the vector mode of linear metric perturbations around spatially flat 
Friedmann-Robertson- Walker (FRW) spacetimes. The procedure we use is analogous to 
that for scalar modes |20] , although specifics certainly change when vector modes are con- 
sidered. The general form of the perturbed metric around the FRW background containing 
only the vector mode is given by 

fl'OO = , 90a = a^Sa , Qab = a? [5ab + Fa,b + Fb^a] ■ (1) 

The perturbation fields and S*" satisfy F"',a= and S"',a= to separate them from 
scalar gradient terms. In other words these divergence free fields describe the vorticity of 
metric perturbations. An index refers to conformal time t], while a,b, . . . refer to co- 
moving spatial coordinates. In a canonical formulation, the distance element is expressed 
in terms of the spatial metric qab, the lapse function N and the shift vector A^", related to 
the spacetime metric through 

^00 = -A^' + qabN'^N'' , goa = qabN^ , Qab = Qab . (2) 

By comparing expression ([T]) for a perturbed spacetime metric to the relation ([2]), one 
can see that in canonical formulations the vector perturbation is generated through the 
perturbations of the shift vector A^"^ and off-diagonal components of the spatial metric 
qab- In particular, the lapse function, being scalar, does not contribute to vector mode 
dynamics. 

2.1 Background 

In canonical gravity, the spatial metric qab plays the role of a configuration variable with 
momenta related to extrinsic curvature 

Kab=^{qab-2D^aNk)) (3) 

in terms of a spatial covariant derivative Da- The lapse function A^ and shift vector A^" 
do not have momenta and are not dynamical since they do not appear as time derivatives 
in the action. They rather play the role of multipliers to constraints which will be written 
explicitly below for the vector mode. 
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However, in view of including quantum gravity effects we do not use metric variables 
but connection variables which follow after a canonical transformation and provide the 
basis for a loop quantization of gravity [21]. In these Ashtekar variables [151 IE] one uses 
a densitized triad Ef instead of the spatial metric, defined as 

E::=\det{e:M (4) 

in terms of the co-triad and its inverse e" which in turn are defined via the spatial 
metric as qab ='■ e.\e\. The canonically conjugate variable to the densitized triad is the 
Ashtekar connection A\:=T\ + iK^, where is extrinsic curvature and 7 > is the so- 
called Barbero-Immirzi parameter [161 122] • The spin connection by definition satisfies 
Da^i = which can be solved as 

K = -^^'e'>id[A^ + lele[d[A]) (5) 

As mentioned, we perturb basic variables around a spatially-flat FRW background. Fur- 
ther, we denote background variables with a bar as in 

Et = p6t , Tl = , K = k5: , N = VP , iV« = 0, (6) 

where p = a? and the spatial metric is qab = o?5ab- (In general, p as a triad rather than 
metric component could be negative, which however we can safely ignore here.) The choice 
of iV = a leads to conformal time. One may emphasize here that for a spatially flat FRW 
spacetime the background spin connection f ^ vanishes, thus A\ = = •ykSl- Moreover, 
k = da/drj = a/a in conformal time. 

2.2 Perturbed canonical variables 

The perturbed densitized triad and Ashtekar connection around a spatially fiat FRW back- 
ground are given by 

E^ = p6t + 6E^ , A', = Tl + ^K: = ^k6: + i6Ti + ^6K) , (7) 

where p and jk are the background densitized triad and Ashtekar connection. At the linear 
level, the general solution for the co-triad corresponding to qab = a^i^ab + Fa,b +Eb,a ) as in 
equation ([T|) is 

el = a[5: + ic,Fa: + C2F\a)] , (8) 

where C1 + C2 = 1. Specific values of Ci and C2 are part of the gauge choice of the triad as a 
set of three vector fields in arbitrary rotation. Using definition (jlj) for the densitized triad 
one can easily compute the expression of the perturbed densitized triad for vector mode 

<5Ef = -p(ciF„" + C2F",0 , (9) 

where we have used the divergence free property, i.e. SISE"- = 0, for vector mode pertur- 
bations; thus no linear term results from the determinant used in defining the densitized 
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triad. By comparing the perturbed spacetime metric from expression ([T]) with the relation 
([2]), one can read off the expression for the perturbed lapse function and shift vector 

6N = , dN"" = 3" . (10) 

Using the spin connection and the general expression of the perturbed densitized triad 
in ([7]), the linearized spin connection is given by 

STi = -JHa^deSE^ . (11) 



2.3 Canonical structure of linearized vector modes 

In a canonical formulation, the Einstein-Hilbert action can be written equivalently (up to 
boundary terms) using the Ashtekar connection and densitized triad as 



5*^^ = I dt I d'x 



SnG-f 



E^CtK - + N^'Ca + NC] 



'12) 



where A*, A^" and A^ are Lagrange multipliers of the Gauss, diffeomorphism and Hamilto- 
nian constraints, explicit expressions of which are written below. Before decomposing the 
symplectic structure according to ([7]), we introduce a cell to render the homogeneous mode 
well-defined. Integrating the first term of fll2p only over a finite box of coordinate volume 
Vq with perturbed basic variables of the form ([7]), we obtain the symplectic structure given 
by the Poisson brackets of the background and perturbed variables as 

{^,P} = ^ , {5Kl{x),8E]{y)} = MG5\x,y)5',8]. (13) 

In deriving these relation, we have used the properties that for vector perturbations 
Sl^SEf = 5^6Kl^ = 0. This provides separate canonical structures for the background 
and perturbations, but these variables will be coupled dynamically. In particular, the ho- 
mogeneous background is dynamical and would receive back-reaction effects at higher than 
linear orders. 



2.3.1 Gauss constraint 

In triad variables, a Gauss constraint appears which generates internal gauge rotations of 
phase space functions because triads whose legs are rotated at a fixed point correspond to 
the same spatial metric. This constraint is given by 

G(A) := ^ d'xA'G, = ^ d'xA' [daE'i + e,/A^„E«) . (14) 

Using the perturbed form of basic variables ([7j) , it can be reduced to 

G(A) = ^ / d^xA* (e., + e.J'UEl) . (15) 
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Since we are working with a background = pSf whose gauge freedom is fixed, the 
multipher A* is already of first order. To derive the expression (fTSl) . we have used the 
definition of the spin connection in terms of the densitized triad Ef, which reduces to 
( ITTl) for the hnearized equations. Internal gauge rotations of phase space functions are 
parametrized by the Lagrange multiplier A* through 6\f = {f,G[A]}. In particular, the 
internal gauge rotations of perturbed basic variables are 

6a{6K:) := {6KI G(A)} = kA'eiJ , 6a{SE^) := {6E^, G(A)} = -pA'e/^" • (16) 

Clearly, the perturbed variables themselves are not invariant under internal gauge rotations 
in spite of the fixed background. However, one may notice already that the symmetrized 
perturbed variables are in fact invariant under internal gauge rotations. Thus, the physical 
quantities must depend only on the symmetrized form of the perturbed basic variables. 
Then, the constants ci and C2 in ([9]) are irrelevant since 

2.3.2 Diffeomorphism constraint 

The diffeomorphism constraint generates gauge transformations corresponding to spatial 
coordinate transformations of phase space functions. Its general contribution from gravi- 
tational variables is given by 

Dam := ^ d'xN^C, = ^ ^ d'xN^ [F^,E^ - A^G,] (17) 

(where the subscript "G" stands for "gravity" to separate the term from the matter contri- 
bution. In general, a matter field would also contribute a term to the diffeomorphism 
constraint, which we leave unspecified here and relate later to the stress-energy tensor.) 
The second term vanishes by virtue of the Gauss constraint, but is necessary to generate 
diffeomorphisms in the form of Lie derivatives of phase space functions along the shift 
vector. Using the expression of the perturbed basic variables ([7]), one can reduce the 
diffeomorphism constraint to 

= ^ / d'x6N^ [-p{d,6K^) - md.SEt)] . (18) 

Here, we have kept up to quadratic terms in the perturbations, noting N"" = 0. 

2.3.3 Hamiltonian constraint 

In a canonical formulation, the Hamiltonian constraint generates 'time evolution' of the 
spatial manifold for phase space functions satisfying the equations of motion. Its gravita- 
tional contribution in Ashtekar variables is 
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Using the general perturbed forms of basic variables and the expression of curvature 
daAl — dbAl^ + e'^jkA^A'^, one can simphfy (fT9l) . Up to quadratic terms it is given by 



Hg[N] 



16nG 



2p 



2k 



(20) 



with 6N = for vector modes. We may mention here that 7 dependent terms drop out of 
the Hamiltonian constraint once the spin connection, the generic form iQ of the densitized 
triad and the extrinsic curvature are used. Again, there is a matter contribution H^n left 
unspecified here. 



2.4 Gauge transformations and gauge invariant variables 

General gauge transformations are determined by a choice of and A^", which gives rise to 
all space-time coordinate transformations. For the vector mode, only the choice A^" = 
with C,^,a = is relevant since the remaining functions N or Na = C,,a would affect only the 
scalar mode. Thus, transformations of interest here are generated only by diffeomorphism 
constraint in a form parametrized by the shift vector as S^f = {/, Dg[A^" = ^°]}- We 
need to consider this only to linear order, using the linearized constraints. The resulting 
transformation for basic variables will then be used to find gauge invariant combinations. 
Alternatively, a canonical formulation allows one to compute gauge invariant observables 
first, and then linearize [23t [2^. 

Under such a gauge transformation, the perturbed densitized triad and extrinsic cur- 
vature transform as 

S^iSE^) = -pd,C , S^i5K^,) = kdae . (21) 

Relating 6E^ to 6F°- and ^A'* to 55*" based on leads to gauge transformations for the 
vector mode functions and S"". 

S^F"" = C , Si:S^ = C , k^^" = kiS" - F") = (22) 
introducing a"" := S"" — F"" as a. gauge invariant variable for the gravitational vector mode. 

2.5 Linearized equations of motion 

General equations of motion are written canonically as 

/ = {/,^} (23) 

for any phase space function / using the total Hamiltonian Ti. For gravity, Ti = Hg[N] + 
Dg[N'^] + Hm[N] + Dm[N'^] with the gravitational contributions Hg and Dg and matter 
contributions Hm and Dm to the Hamiltonian and diffeomorphism constraints. Equations 
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of motion refer to coordinate time, with derivatives being indicated by the dot. The form 
of the lapse function N specifies which time coordinate is used; here, N = a = ^J^p imphes 
conformal time r/. The general form (|23l) also applies to equations of motion for momenta 
of the multipliers, such as P/v = = {P/v,7Y} = —5T-i/5N. This must be zero because the 
momentum P/v is zero for an action ( JT2l) not depending on the time derivative of A^. In 
this way, equations of motion for the momenta of lapse and shift give rise to constraints. 
Hamilton's equation of motion for the perturbed densitized triad 



leads to the expression ([3]) 
symmetrized form which is 



5E^ = H} (24) 
of extrinsic curvature, but linearized. We will only need it in 



or 



using the gauge- invariant variable a". Variation with respect to the shift vector 5A^", 
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(25) 
(26) 

(27) 



when expressed in terms of symmetrized vector perturbations implies 



SttG 



f28) 



Using the relation between the perturbed stress-energy tensor and a variation of the matter 
diffeomorphism constraint Dm with respect to the shift vector as derived in [20] . 



]Vp3/2 



(29) 



one can express equation fl28|) in standard form [25 

1 



2a2 



VV, = 87rGip + P)Va 



(30) 



for a vector mode equation. Here p and P are energy density and pressure of the background 
matter field. 

The second vector mode equation comes from Hamilton's equation of motion for per- 
turbed extrinsic curvature, which using (12^ becomes 



1 d 
2d^ 



a) 



(31) 
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The perturbed spatial stress tensor for the vector mode in terms of the matter Hamiltonian 
10] is 



P 



1 OH fSE^5^, 5f 



3Vo dp 



P 



+ 



(32) 



By using the expression of the perturbed stress-energy tensor in terms of anisotropic stress 
TTa as TTa,* + 7r*,a ) we cau express equation fl3T|) in standard form [25] 



(33) 



for the second vector mode equation. 

For vector modes, variation with respect to the lapse function does not give new field 
equations but would rather contribute back-reaction terms to the background evolution. 
The two equations (|28l) and (1311) thus provide the complete dynamics. 



3 Quantum dynamics: Inverse triad corrections 

We have completed the derivation of vector mode equations in Hamiltonian cosmology 
based on Ashtekar's formulation of general relativity. This naturally agrees with results 
of [25]. As mentioned, we are interested in applying this formulation to study possible 
canonical quantum gravity effects. As an explicit example, we now consider quantum 
corrections coming from terms containing inverse densitized triad components. In Sec. H] 
we provide formulas for a second major quantum correction that one expects from the 
use of holonomies, rather than direct connection components, as basic operators of the 
quantum theory. This has the effect of adding terms of higher order in extrinsic curvature 
components, and thus higher powers of the first time derivative of the metric, to the 
Hamiltonian. In addition to those two corrections, there are higher derivative corrections 
implied by genuine quantum effects [261 122] • AH this combines to effective constraints or 
effective equations of motion for the system. As for the diffeomorphism constraint, we 
assume that it receives no quantum corrections because it is quantized directly through its 
phase space transformations |28j . 



3.1 Quantum corrected Hamiltonian constraint 

While homogeneous quantum cosmology using loop quantum gravity techniques is rather 
well understood [H], a systematic derivation of quantum corrections to classical dynam- 
ics which includes inhomogeneity is not yet available. But typical effects are known and 
provide valuable indications for implications of quantum effects. In loop quantum gravity, 

the appearance of inverse powers of the densitized triad as in , initially leads to 

\/\det E\ 

difficulties since flux operators quantizing the densitized triad have discrete spectra con- 
taining zero as an eigenvalue [29], [30]. These difficulties can be overcome in a way exploiting 



9 



background independence, and giving rise to well-defined operators [311 [32]. However, for 
small values of the densitized triad, where the classical expression would diverge, there are 
deviations from the classical behavior which imply quantum corrections. For the homo- 
geneous case, explicit calculations show that the classical term is multiplied by a factor 
a [33] making the whole expression finite. For large fluxes, the leading terms in explicit 
expressions of a are of the form 

a(p) = l + c(^^j , (34) 

where n and c are positive deflnite numbers which correspond to a given inverse triad oper- 
ator. However, they are not completely flxed since triad operators themselves are subject 
to quantization ambiguities [3H [35] . One important motivation to study inhomogeneous 
models is that, compared to homogeneous models, their dynamics gives rise to much tighter 
consistency conditions which could constrain such parameters. This is indeed borne out, 
as we will see later. While general derivations of a for inhomogeneous conflgurations with- 
out assumptions on the metric such as symmetries or speciflc modes is complicated (see 
e.g. [3B]) and does not yet provide many insights, a can be computed in certain pertur- 
bative regimes. It has been studied recently for scalar mode perturbations [20j, showing a 
similar behavior as in the isotropic case. Here we consider such quantum corrections for 
vector mode perturbations, although this is a case where no explicit expression is available 
yet. As we will see, consistency itself restricts the form of a beyond what can currently 
be computed directly from operators. Thus, our perturbative treatment at the effective 
level provides conjectures to probe the overall consistency of the theory by comparing with 
results for the underlying operators. 

This procedure leads us to an ansatz for the quantum corrected Hamiltonian constraint 



2k 



(35) 



where a{p, SE°-) is the correction function, now also depends on triad perturbations. It 
is important to emphasize here that the correction a coming from the quantized inverse 
densitized triad in general could be tensorial in nature. However, later we will see that 
the leading effect on perturbation dynamics comes from the background corrections, i.e. 
from a = a{p,0), as well as derivatives of a evaluated at the background conflguration. 
(Note that the only background variable determining the geometry is p, in which the 
corrections are expressed. This is sometimes seen as problematic since the scale factor 
a = a/P can be rescaled arbitrarily in a flat isotropic model. However, the dependence of a 
function a in an inhomogeneous Hamiltonian constraint is through elementary fluxes whose 
values are determined by an underlying inhomogeneous state. The scale of corrections, 
too, is determined by the underlying state, resolving any apparent contradiction between 
the appearance of such a scale and the rescaling freedom of a flat, precisely isotropic 
background.) 
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3.2 Constraint algebra 



Ac^ = 3Nk'y/p 



In a canonical formulation of general relativity, classical constraints form a first class Pois- 
son algebra, i.e. {Ci,Cj} = flj{A, E)Ck, whose coefficients can in general be structure 
functions. It ensures that the transformations generated by the constraints are gauge and 
preserve the constraint surface. In other words, the evolution of phase space functions pre- 
serves the physical solution surface. To study quantum gravity effects, we have introduced 
a quantum correction function a{p, SE^) which depends on phase space variables. Natu- 
rally, having a new expression for the Hamiltonian constraint, there could be an anomaly 
term of quantum origin in the constraint algebra. Here, while {Hq[N\, Hq[N'W is trivial 
in the absence of lapse perturbations for the vector mode, a non-trivial anomaly in the al- 
gebra could occur in the Poisson bracket between Hq[N] and DqIN""]. This bracket turns 
out to be 

{H^iNlDcm) = d\p{d,6N'^)A.^ , (36) 

where 

We see that the anomaly term contains derivative of a with respect to both p and SE^. 
However, as mentioned before the functional form of the correction function a(p, 5E^) in 
terms of 5E^ is not known while the p-dependence can be taken to be of scalar mode form, 
i.e. with parameters c and n fixed once an inverse triad operator is chosen. To have a 
consistent set of evolution equations for the vector mode we require the anomaly term to 
vanish i.e. Ac = 0. This in turn puts restrictions on the linearized functional form of a as 
a function of 5E°;: 

da 1 da frj^dckcj\ 

WW) ^ ~^p9P ^ ■ 

Since a is in principle computable in the full theory, this provides important consistency 
checks for loop quantum gravity. At present, only the dependence a{p) as well as derivatives 
of a along diagonal components of the spatial metric are known . Anomaly cancellation 
will then lead us to conjecture a form of derivatives Sa/S{SEf) along off-diagonal compo- 
nents of the metric which one can later compare with direct calculations once they become 
available. 



3.3 Effective gauge invariant perturbation and its linearized equa- 
tion of motion 

Extrinsic curvature is derived using one of Hamilton's equations of motion. Thus, one 
expects the expression of extrinsic curvature as it follows from an equation of motion to 
change due to the quantum corrections. This incorporates an effect of quantum geometry 
which changes the usual differential geometric relation of extrinsic curvature corresponding 
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to changes of the spatial metric between different shces. One can easily compute the 
corrected expression for the perturbed extrinsic curvature 

5i^(a^^ = -^(^a/ + ) + + ) (38) 

where a appears only in the first term expressed through the classical gauge invariant 
quantity = S'' - F'^. 

The gauge transformation of the vector perturbation remains unchanged because 
the diffeomorphism constraint retains its classical form. Even though there is a quan- 
tum correction to extrinsic curvature one can easily see that a" = S*" — is still the 
gauge invariant variable. Moreover with the diffeomorphism constraint being unaffected, 
Eq. ( l28i) remains unchanged for the quantum background dynamics. Now using the second 
Hamilton's equation, one obtains an equation of motion for extrinsic curvature as follows 

+ A(^^ = STxGp5T^''^^\) . (39) 

To derive Eq. (1391) , we have used the corrected expression ( |38l) of extrinsic curvature as well 
as an analogous expression for the background extrinsic curvature, ok = a/ a. One may 
note here that equation fl39|) explicitly contain the anomaly term. Thus, the requirement of 
an anomaly free constraint algebra leads to the corrected equation fl39|) explicitly in terms 
of the gauge invariant variable. The presence of an anomaly, on the other hand, would 
make it impossible to express the equations of motion solely in terms of gauge invariant 
variables. Since consistency of the constraint algebra requires us to set = 0, closed 
equations for the gauge invariant perturbations follow. Nevertheless, non-trivial quantum 
corrections remain through a. 

Note that the only correction function a in the Hamiltonian appears in a form multiply- 
ing the lapse function N . The correction could thus be reduced to a simple change of the 
lapse function and thus the time gauge. Still, the corrections are non-trivial as illustrated 
by the equations of motion shown here. The choice and interpretation of time is based 
on the line element since this determines the measurement process of co-moving geodesic 
observers. This behavior is not changed by the appearance of a correction function a in 
the Hamiltonian constraint even if it always appears in combination with N . Even with a 
corrected Hamiltonian constraint we are still referring to the same form of conformal time, 
but fields evolve differently as given by the corrected constraint. Then, also observable 
implications of the quantum corrections are possible. 

4 Quantum dynamics: Holonomy corrections 

In addition to corrections in coefficients of the constraint due to inverse powers of densi- 
tized triad components, there are corrections which resemble higher curvature terms in an 
effective action. While these corrections would be dominant in purely isotropic models by 
virtue of the large matter energy density in a macroscopic universe, they are sub-dominant 



1 

a 



1 d 
'2d^ 



+ CT\a )-k{a- a'p){aa' + a\a , 
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in inhomogeneous situations [37]. Moreover, as we will see they do not provide much of 
a structural change to the equations. The Hamiltonian constraint operator is formulated 
in terms of holonomies rather than connection or extrinsic curvature components. Since 
these objects are non-linear as well as (spatially) non-local in connection components, they 
provide higher order and higher spatial derivative terms. Higher time derivatives, as they 
would also be provided by higher curvature terms, do not arise in this way but rather 
through the coupling of fluctuations and higher moments of a quantum state to the expec- 
tation values [261 122]. The full effective constraint including all these terms has not been 
derived yet. In this section we therefore focus on an analysis of higher order terms only. 

For an isotropic model sourced by a massless, free scalar field such higher order terms 
turn out to be the only corrections, and there are no higher time derivatives [38]. The exact 
effective Hamiltonian can then be obtained by simply replacing the background Ashtekar 
connection 7A; by sin fi'-fk, as it was also seen in numerical studies [39]. The parameter 
/i depends on the quantization scheme and may be a function of p. Just as with the 
parameter c in fl34p . we will see that the freedom is constrained by anomaly cancellation. 
To study the effects of the background dynamics on inhomogeneous perturbations, we 
similarly substitute the appearance of k in the classical Hamiltonian by a general form 
{mfl)~^ sin mjl'-fk where m is an integer. (This parameter is kept free because different 
factors of sines and cosines combine from the full constraint to result in this term. It can 
be constrained by looking at detailed properties of the underlying operator, but also by 
consistency requirements as we will see shortly.) With this prescription, one can write 
down expression for the corrected Hamiltonian constraint 



2pi 



(40) 



Here we have required that the effective Hamiltonian ( HOl) has a homogeneous limit in 
agreement with what has been used in isotropic models. This fixes the parameter m to 
equal one in the first two terms. The parameter for the last term as chosen here is the one 
which leads to an anomaly-free constraint algebra. 

Although we write explicit sines in this expression, and thus arbitrarily high powers 
of curvature components, it is to be understood only as a short form to write the leading 
order corrections. Higher orders are supplemented by further, yet to be computed higher 
curvature quantum corrections. The expressions are thus reliable only when the argument 
of the sines is small, which excludes the bounce phase itself. (Such sine corrections can 
be used throughout the bounce phase only for exactly isotropic models sourced by a free, 
massless scalar [3HI SO] , but not in the presence of a matter perturbation [H] or anisotropics 
and inhomogeneities.) 
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4.1 Constraint algebra 



Again, a non-trivial anomaly in the algebra can occur between the Poisson bracket between 
H'^[N] and Dg[N''] 



{H^[NIDg[N'^]} 



50 



sin 2^'-)k 

'^m ) 



where 



N 

Vp 



_ d f sin fijk 
^ dp 



/i7 



Dg[N-] + 



f sin /i7A; 



SvrG 



/6E: 



V P 



(41) 



(42) 



One may easily check here that p ~ 1/v^ leads to an anomaly free algebra up to order ¥. 
This is in accordance with the result of arguments put forward recently in purely isotropic 
models jl2]. From an inhomogeneous perspective, the behavior /2 ~ I/a/p reflects the fact 
that the fundamental Hamiltonian creates new vertices when acting on a graph state such 
that the number of vertices increases linearly with volume US]- This suggests a tight 
relation between anomaly freedom at the effective level and properties such as the creation 
of new vertices by a fundamental Hamiltonian constraint. 



4.2 Effective linearized equation 

Hamilton's equations governing the background dynamics are given by 



and 



Vp 



p 



sin fi'yk^ ^ 



_ f sin 2fi'~fk\ 



2p( 



V 2^7 ; 



+ Pt- 



— 2 

d I sin p'-fk 



+ 



SttG fdH, 



SVo V dp 



(43) 



(44) 



Since the diffeomorphism constraint is assumed to remain unaffected, the gauge transfor- 
mation of vector functions remains unchanged. However, a correction to the extrinsic 
curvature expression now leads to a new expression for the gauge invariant variable 



a^ = S'^-F^ + k(l-'^^] 



(45) 



The expression of perturbed extrinsic curvature in terms of the gauge invariable variable 
remains unchanged, though. 



(46) 



Now using again Hamilton's equation, one obtains an equation of motion for extrinsic 
curvature as follows 



1 d 
2d^ 



(47) 



As before, an anomaly-free constraint algebra, requiring A\ = 0, leads to a quantum 
corrected equation entirely in terms of the gauge invariant variable. 
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5 Rate of change of vector perturbations 



We now have equations that govern the dynamics of vector mode perturbations including 
quantum corrections. For simphcity, we consider the situation where anisotropic stress is 
absent, which is the context of [25| . 

5.1 Classical dynamics 

Let us recall the key feature of classical vector mode perturbations. In the absence of 
anisotropic stress, the right hand side of equation (pTl) vanishes. The classical background 
extrinsic curvature is related to the time derivative of the scale factor a as k = a/a. It 
is also convenient to decompose the amplitude of perturbations in terms of their Fourier 
modes cx^. With these simplifications, equation ( pTll leads to the rate of change 



for Fourier modes. Thus, any vector mode grows as al ~ in a contracting phase, and 
correspondingly decays in an expanding phase. This is independent of the background 
matter content except for the assumed absence of anisotropic stress. Thus, a quantum 
correction of the background matter sector alone would not help in taming the growth 
of vector perturbations. Intuitively, one expects that it is required to have a modified 
gravity sector in order to have any modification of the growth rate. In the next sub-section 
we consider the dynamics in the presence of quantum corrections to the gravitational 
Hamiltonian. 

5.2 Quantum corrections: Inverse triad 

As in the classical case, we consider the situation where anisotropic stress is absent. Then 
using the equation of motion ak = a/ a for the background, one can write Eq. ( l39l) for 
Fourier modes al of vector perturbations in the form 



For a = 1 quantum corrections are switched off and we obtain the classical result. Now 
using the generic form of a as in flMl) with an approach to one from above (i.e. c > 0), it 
is easy to see that 



Thus, the decay rate of vector mode perturbations is slightly higher compared to that of 
a classical scenario. In other words, in the contracting phase, the correction coming from 
inverse powers of the densitized triad in background dynamics causes vector perturbations 
to grow even faster, though only slightly, than in the classical situation. The quantum 




(49) 




(50) 
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correction depends inversely on volume, i.e. it becomes stronger in the smaller volume 
regime. Using non-perturbative corrections in a for small densitized triads, a decrease of 
the rate is indicated since a falls below one in this regime and has a' > 0. However, for such 
small scales the perturbation theory of inhomogeneities is less reliable and a suppression 
of the decay rate on very small scales can, at present, at best be taken as an indication. 

5.3 Quantum corrections: Holonomies 

As before, we consider the situation where anisotropic stress is absent. Now, the corrected 
equation for Fourier modes a], of vector perturbations, after dividing (H7j) by (H3l) . is 



Here, quantum corrections disappear for /i — 0. Again, the right hand side is less than 
—2 and thus vector modes grow more strongly in a contracting phase. If the behavior 
is extrapolated to the bounce phase, the growth rate in a contracting universe becomes 
even larger and would diverge at the bounce where cos{fl'jk) is zero. This indicates a 
breakdown of the perturbation scheme and the need to include higher order terms as the 
bounce is approached. It should also be emphasized here that /i ~ 1/ ^/p leads to anomaly 
free Poisson algebra only up to order k^. So the anomaly term becomes significant near 
bounce phase that makes the analysis less reliable there. 

6 Discussions 

In the absence of anisotropic stress, gauge-invariant vector perturbations classically grow 
as in the contracting phase. Such a growth of vector perturbations indicates a pos- 
sible violation of homogeneity assumptions in smaller volume regime, indicated by the 
breakdown of classical perturbation theory. Thus, the growth of vector perturbations may 
pose significant problems in particular for bouncing cosmologies which are invariably as- 
sociated with a contracting phase but often have been derived only under the assumption 
of homogeneity. In these models, conclusions regarding bounces are drawn based on the 
homogeneity assumption. Naturally, a growth of vector perturbations can question the 
robustness of such bounce scenarios by questioning the validity of the homogeneity as- 
sumption itself at smaller volume, given that bounces are typically more difficult to realize 
when inhomogeneities are taken into account; see [SIEH]- Thus, it is an important issue 
to study the dynamics of vector mode perturbations in the cosmo logical context. 

In this paper, we have presented a systematic derivation of gauge-invariant vector per- 
turbation equations to linear order in Hamiltonian cosmology based on Ashtekar variables. 
We have only considered a spatially flat Friedmann-Robertson- Walker background as this 
is the case of most interest in cosmology. Hamiltonians and equations of motion are techni- 
cally more complicated in the presence of spatial curvature and are still being worked out. 
Quantum corrections are, however, analogous in spatially curved backgrounds and we do 




(51) 
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not expect our results to change significantly in those cases. Specifically, we have studied 
the effects of two particular types of quantum correction, inverse triad and holonomy cor- 
rections, on the dynamics of vector perturbations in large volume regimes. For each type, 
we have shown that in a contracting phase the growth rate of vector mode perturbations is 
slightly stronger compared to the classical situation due to quantum effects. This quantum 
correction is small as one expects. Although there are quantization ambiguities, the sign 
of corrections to the growth rate seems robust. A reduction of the growth rate is indicated 
only in regimes of non-perturbative corrections of inverse powers. For such a reduction to 
be realized, even if it remains true under a more careful perturbation analysis, one would 
have to enter the deep Planck regime. Moreover, if one starts with a large classical universe 
as initial configuration, such non-perturbative quantum effects will become relevant only 
after long evolution times. In general, due to the growth of the vector mode one will even- 
tually have to use higher than linear orders in pcrturbative inhomogeneities which we have 
not included in this paper. This by itself may well change some of the conclusions about 
the bounce phase independently of quantum effects in the evolution of inhomogeneities. 

Another important issue touched in this paper is that of potential anomalies in the 
quantum constraint algebra. We started with a general but unspecified form of a quantum 
correction function a or higher order terms, including inhomogeneity. The functional 
form of q; as a function of the background variable p is known. However, its functional 
dependence on the perturbed densitized triad bE^ (which is purely off-diagonal for vector 
modes) is unknown due to the lack of a systematic derivation of such corrections from the 
full theory. As we have observed, requiring anomaly cancellation in the modified constraint 
algebra restricts the functional dependence of quantum correction functions such as a on 
off-diagonal triad components. It would be interesting to see whether such a restriction 
is satisfied by a systematically derived quantum correction function from the full theory. 
As a key result, we have observed the possibility of non-trivial quantum corrections while 
preserving anomaly freedom. The classical constraint algebra for vector modes is rather 
trivial, but is much more restrictive for scalar modes for which the calculations here show 
the guiding principle. 
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